Excitation of inertial modes in a closed grid turbulence experiment under rotation 
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We report an experimental study of the decay of grid-generated turbulence in a confined geometry 
submitted to a global rotation. Turbulence is generated by rapidly towing a grid in a parallelepipedic 
water tank. The velocity fields of a large number of independent decays are measured in a vertical 
plane parallel to the rotation axis using a corotating Particle Image Velocimetry system. We first 
show that, when a "simple" grid is used, a significant amount of the kinetic energy (typically 50%) 
is stored in a reproducible flow composed of resonant inertial modes. The spatial structure of those 
inertial modes, extracted by band-pass filtering, is found compatible with the numerical results of 
Maas [Fluid Dyn. Res. 33, 373 (2003)]. The possible coupling between these modes and turbulence 
suggests that turbulence cannot be considered as freely decaying in this configuration. Finally, 
we demonstrate that these inertial modes may be significantly reduced (down to 15% of the total 
energy) by adding a set of inner tanks attached to the grid. This suggests that it is possible to 
produce an effectively freely decaying rotating turbulence in a confined geometry. 



I. INTRODUCTION 

Translating a grid in a closed volume of fluid is a stan- 
dard way to generate an approximately homogeneous and 
isotropic turbulence in order to investigate its temporal 
decay. Although the level of homogeneity and isotropy 
of the turbulence generated in this way is not as good as 
in the more conventional configuration of a fixed grid 
in an open wind tunnelj^l this closed flow configura- 
tion has proved to be very useful when a compact sys- 
tem is needed, and in particular when experiments are 
performed in a rotating frameP^ Apart from early ex- 
periments performed in wind tunnels with a rotating 
sectionj^ all rotating grid-generated turbulence exper- 
iments since then are based on oscillated or translated 
grids in closed containers. 

Although grid turbulence in open and closed geome- 
tries share similar properties, large scale recirculations 
are more likely produced in the closed configuration.^ 
These reproducible flows, which will be called mean flows 
in the sequel, can be easily extracted thanks to ensemble 
average over several independent realizations. They may 
actually be triggered from inhomogcncitics in the wake 
of the grid, or from spontaneous symmetry breaking in 
a perfectly symmetric geometry. In the case of decaying 
turbulence, although the kinetic energy of the mean (i.e., 
ensemble-averaged) flow may be negligible compared to 
that of turbulence soon after the grid translation, it may 
persist over large times and possibly influence the final 
period of decay of turbulence. 

In the presence of background rotation, grid-generated 
turbulence may excite inertial waves. These waves are 
anisotropic transverse dispersive waves, which propagate 
through a rotating fluid because of the restoring nature 
of the Coriolis force Their pulsation a lies in the 
range [0, 2S7], where 17 is the rotation rate, indicating that 



these waves can be excited when the characteristic time 
of the considered flow is of order of the rotation period, 
i.e. when the Rossby number is of order unity. These 
inertial waves may originate either from a mean flow, in 
which case they are reproducible and can be detected in 
the ensemble-average, or from turbulence itself, in which 
case they are non reproducible and cancel out by phase 
mixing in the ensemble-average. 

In grid-generated turbulence experiments performed in 
a closed container, still in the presence of rotation, re- 
producible osc illatio ns induced by inertial modes have 
been reported ! 4 * 6 * 13 ! A first mechanism responsible for 
these modes is the resonance of multiply reflected inertial 
waves over the containerP^E^ Another possibility is that 
they originate from direct conversion of the energy of the 
large scale recirculations initiated by the grid translation, 
which were also present in the non-rotating case, when 
the characteristic time of these recirculations decreases 
down to the order of the rotation period. 

Inertial modes can be found in general in a container 
when the walls are either normal or parallel to the rota- 
tion axis, whereas wave attractors are found when sloping 
walls are present.^ Only the former are considered in the 
following. Their resonant frequencies can be derived an- 
alytically only in some specific geometries, such as the 
so-called Kelvin modes in a cylinder rotating about its 
symmetry axisP^ In the case of a parallelepipedic con- 
tainer, the frequencies and the spatial structure of the 
eigenmodes have been characterized in detail by Maas^l 
for an inviscid fluid from numerical simulations. 

Inertial modes excited in grid-generated turbulence 
have been first observed in a parallelepipedic channel 
with a free surface by DalzielP who pointed their poten- 
tial influence on the decay of the turbulent component of 
the flow. Inertial oscillations are also clearly visible in the 
experiments of Morize and Moisy-- (with a rigid upper 
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surface) and Moisy et oiP^ (with a free surface). They 
have been characterized by Bewley et alW in a set of two 
experiments, in which a grid is rapidly drawn in liquid 
helium or nitrogen in cylindrical and squared geometries. 
These authors found good agreement between the mea- 
sured frequencies and the numerical results of Maas^l 
for various aspect ratios. They conclude that translat- 
ing a grid in a closed rotating container does not allow 
to generate freely decaying turbulence, because a signif- 
icant amount of the energy goes into inertial modes. In 
all these experiments, these inertial modes are detected 
in the ensemble-averaged flow, suggesting that they are 
more likely produced by reproducible flow features, than 
from turbulence. 

The aim of the present paper is to investigate in more 
details the structure of the inertial modes excited when 
towing a grid in a rotating square container, and to ex- 
plore to what extent those modes may be redu ced. An 
experimental setup similar to that of Morize et alW^ has 
been mounted on a new rotating platform, allowing for 
Particle Image Velocimetry measurements in a vertical 
plane. We demonstrate that, by attaching a set of in- 
ner sidewalls to the grid, the amount of energy stored in 
the inertial modes may be drastically reduced. A simi- 
lar configuration has been first shown by Staplehurst et 
aZP to reduce significantly the large-scale recirculations 
in the absence of rotation. We demonstrate here how a 
configuration inspired from their work inhibits the exci- 
tation of inertial modes, by an enhanced conversion of 
the mean flow energy into turbulence. 



II. EXPERIMENTAL SETUP 
A. Tank and rotating platform 




FIG. 1. (Color online) Schematic view of the experimental 
setup, in the "simple grid" configuration. The 60 x 60 x 60 cm 3 
tank is filled with 52 cm of water and is rotating at an angular 
velocity of fi = 0.84 rad s _1 . A square grid of 40 mm mesh is 
towed by four shafts from the bottom to the top at constant 
velocity V g = 0.70 m s . PIV measurements in a vertical 
plane (x, z) are achieved in the rotating frame, based on a 
laser sheet illuminating the vertical plane (a;, z) and a camera 
aiming normally at it. 



brushless motor ensuring the vertical translation of the 
grid, with acceleration and deceleration phases less than 
18% of the total translation time. The Reynolds number- 
based on the grid mesh is Re g — V g M/v = 28 000, and 
the Rossby number is Ro g = V g /2QM = 10.4, indicat- 
ing that the flow in the close wake of the grid is fully 
turbulent and weakly affected by the rotation. 



Only the "simple grid" configuration, used in Sec. Ill 
is described here. The modified configuration with addi- 
tional inner walls is detailed in Sec. HVl 

The experimental setup consists in a cubic glass tank, 
of lateral side L — 60 cm, filled with 52 cm of water (see 
Fig. [T]) and mounted on a precision rotating turntable 
of 2 m in diameter ! 12 * 21 ! The angular velocity f2 of the 
turntable is set to 0.84 rad s _1 (8 rpm), with relative 
fluctuations AJ7/f2 less than 5 x 10~ 4 . The rotation of 
the fluid is set long before the experiment (at least half an 
hour) in order for transient spin-up recirculations to be 
damped and therefore to achieve a solid body rotation. 
A horizontal cover is placed at a height of H = 49 cm, 
defining the upper boundary of the flow. 

Turbulence — and inertial modes — is generated by 
rapidly towing a square grid at a constant velocity V g = 
0.70 m s _1 from the bottom to the top of the tank. Dur- 
ing the subsequent decay of turbulence, the grid is kept 
fixed, at a height of 46 cm, slightly below the cover. The 
grid consists in 8 mm thick square bars with a mesh size 
M = 40 mm, and a ratio of solid to total area of 0.44. It 
is rigidly attached at its four corners to a servo-controlled 



B. PIV measurements 

Velocity fields in a vertical plane (x, z) are measured 
using a 2D Particle Image Velocimetry (PIV) system.^ 
The flow is seeded with 10 /um tracer particles, and il- 
luminated by a corotating vertical laser sheet passing 
through the center of the tank and generated by a 140 mJ 
Nd:YAg pulsed laser. The entire 60 x 46 cm 2 flow sec- 
tion is imaged through a transparent side of the tank 
with a double-buffer high-resolution 2048 x 2048 pixels 
camera, corotating with the tank and aiming normally at 
the laser sheet. During the decay of turbulence, 428 im- 
age pairs are acquired at a sampling rate of 2 Hz. Since 
the typical flow velocities decrease with time, the delay 
between the two successive images of a pair is made to 
gradually increase during the acquisition sequence, from 
10 to 68 ms, so that the typical particles displacement 
remains constant, of order of 5 pixels, during all the de- 
cay. PIV computations are then performed over image 
pairs, on 32 x 32 pixels interrogation windows with 50% 
overlap, leading to a spatial resolution of 4.9 mm. 
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FIG. 2. (Color online) Vertical velocity u z (xo, zo,t) at the 
center of the flow for 20 realizations performed at Q = 
0.84 rad s" 1 (in various colors). We superimposed in black 
thick line the ensemble average of these time series. 

C. Reynolds decomposition 

In Fig. [2j we show the time series of the vertical ve- 
locity u z (xq, zo,t) at the center of the flow for 20 inde- 
pendent realizations of the decay. The origin t = is 
defined as the time at which the grid reaches the top of 
the tank. The ensemble-average of those realizations is 
also shown in bold line. This plot clearly illustrates that 
the flow consists in well-defined oscillations excited by 
the grid, of characteristic timescale of about one period 
of rotation, superimposed to non-reproducible turbulent 
fluctuations. Those oscillations actually correspond to 
inertial modes, and their amplitude is clearly of the or- 
der of the turbulence. 

In order to investigate properly the dynamics of the 
inertial modes and of the turbulence, and the possible 
coupling between the two, we introduce the standard 
Reynolds decomposition of the velocity: 

u(x,f) =U(x,t) + u'(x,i). (1) 

Here u(x,i) is the total velocity field, U(x, i) = u(x, t) 
its ensemble average (i.e. the reproducible component 
of the flow) and u'(x, t) its turbulent component. The 
overbar • stands for the ensemble average over several 
independent realizations of the flow. Note that, although 
the ensemble average is often approximated by a tempo- 
ral or a spatial average in most turbulence experiments, 
the use of true ensemble averages here is critical to sepa- 
rate properly the reproducible non-stationary component 
of the flow from the turbulence. 

The Reynolds decomposition ([I]) naturally leads to in- 
troduce three kinetic energies, characterizing the total, 
mean and turbulent flows respectively, and defined as: 

k tot (t) = (u^x,*)), 

kmean{t) — (U ( x j^))) 

fc t „ r6 (t) = ^(x,i)), (2) 




FIG. 3. (Color online) Total (dashed), mean (continuous) 
and turbulent (dashed-dotted) kinetic energies as a function 
of the number of tank rotations Q,t/2n from 40 realizations 
performed at Q — 0.84 rad s _1 . Inset: Ratio of turbulent to 
total kinetic energy a Q, measured at times t n of maximum 
mean energy. 



satisfying k tot (t) = k mean (t) + k tur b{t) ■ Here the brackets 
(•) denote the spatial average over the whole fluid volume. 

In practice, only the two velocity components u x and 
u z can be measured, the measurements being restricted 
to a vertical plane (x, z) located at mid-width of the 
tank. The measured kinetic energies must therefore be 
considered as approximations of the true ones. For a 
statistically homogeneous and isotropic turbulent com- 
ponent, one would simply have k^ r s b = 2k tur b/3. The 
situation is more complicated for the energy of the mean 
flow for two reasons: (i) during one period of a given in- 
ertial mode, the kinetic energy is basically exchanged be- 
tween the measured and the non-measured velocity com- 
ponents; (ii) the ratio between the energy averaged over 
the whole tank and the energy averaged over the mea- 
surement plane (x, z) only depends on the details of the 
spatial structure of each mode, which are not known a 
priori. For those reasons, we do not apply here any cor- 
recting weight to the velocity components when comput- 
ing kinetic energies, and we simply define 

ktot(t) = (u x 2 (x, z, t) + u z 2 (x, z, t)) x>z (3) 

(with {-) x ,z the average over the vertical plane), and simi- 
larly for k mean and k tur b- Note that those modified defini- 
tions still satisfy the relation k to t(t) — k mean (t)+kt U rb(t). 

III. INERTIAL MODES PRODUCED BY THE 
SIMPLE GRID CONFIGURATION 

A. Kinetic energy decay 

In Fig. [3j we present the time evolution of the three 
kinetic energies (pi) computed from 40 decay realizations 
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performed at £1 = 0.84 rad s _1 . Both the energy of the to- 
tal and the mean flow show, superimposed to their overall 
decay, marked oscillations corresponding to the inertial 
modes. On the other hand, the turbulent energy shows a 
monotonous decrease, suggesting a good separation be- 
tween the reproducible and non-reproducible components 
of the flow. At this point, it is important to note that the 
kinetic energy of the mean flow k mean (t) is of the same 
order than the turbulent one kt U rb(t). 

The oscillations of fc tot (i) and k mean (t) are only due 
to the 2-dimensions and 2-components restriction of the 
PIV measurements. Indeed, a monotonous decay should 
be expected for the true energies ^ computed from the 
3 velocity components averaged over the whole fluid vol- 
ume. We expect therefore that the total and mean kinetic 
energies are correctly estimated only at the maximum 
of their oscillations. In order to evaluate the relative 
amount of turbulent and mean energy in the total flow, 
we shall therefore consider only the times t n at which 



l (t) is maximum, and we introduce the ratio: 



a(t n ) 



(tn) 
ktot(tn) 



(4) 



In the inset of Fig. [3j we actually see that turbulence 
represents only 50 ± 10% of the energy in the flow. This 
quite low ratio indicates that the turbulence produced 
with this "simple" grid configuration is evenly distributed 
among the reproducible inertial modes and the "true tur- 
bulence" . It is therefore questionable to consider this tur- 
bulence as freely decaying. Indeed, the possible coupling 
between the turbulence and the mean flow, which cannot 
be investigated at this point, may prevent the turbulence 
to decay freely, as it is continuously fed through energy 
transfer from the inertial modes. 



B. Fourier analysis of the inertial modes 

In order to characterize in more details the structure of 
the flow generated by the grid translation, we perform a 
temporal Fourier analysis of the ensemble-averaged How 
U(x, t) = u(x, t). Assuming that this mean flow is com- 
posed of inertial modes only, U(x, t) can be written as 
follows: 



U(x,t) = 5ft ^ a nms (t) v nms (x) e l 



(5) 



where v„ ms (x) is the (complex) spatial structure of the 
[n, m, s] mode and w„ ms its pulsation, which lies in the 
range [0, 20] allowed for inertial waves. Here, ^t(X) 
stands for the real part of the complex number X. With 
notations similar to that of Maas^ and Bewley et al., 6 
we label here the modes using two integer indices, n and 
77i, and a sign, s — ±. The first index n is the normal- 
ized vertical wavenumber such that the horizontal (resp. 
vertical) velocity component has n (resp. n — 1) nodes in 



the vertical direction. The horizontal structure of a given 
mode of vertical index n is characterized by the second 
index to. Larger values of to essentially correspond to 
finer structures in the horizontal direction, although it 
is not directly related to the number of nodes as for the 
vertical index n. Finally, the sign s refers to the symme- 
try of the mode with respect to the rotation axis: s = + 
for a symmetric mode and s = — for an antisymmetric 
mode. The pulsations uj nms are increasing functions of n 
and, at fixed n and s, decreasing functions of to. In the 
absence of coupling with other modes or with turbulence, 
the amplitude of each mode, a nms (t), is expected to be a 
decreasing function of time because of viscous damping. 

The modes present in the mean flow are identified from 
the temporal Fourier transform of the ensemble-averaged 
velocity field U at each position in the (x, z) measure- 
ment plane, 



U a (x,z) 



TJ(x, z, t) e 



dt. 



(6) 



The temporal bounds T m j„ and T max have been chosen 
equal to 1 and 15 rotation periods, respectively. From 
this, we define the temporal energy spectrum of the mean 
flow, spatially averaged over the (x, z) plane, as: 



E,, 



1 

2^ 



(|U CT (x,z)| 2 



(7) 



In Fig.|4j the spectrum E mean (a) clearly shows a series 
of peaks, whose values are listed in Table [TJ These peak 
frequencies are in excellent agreement with some of the 
numerical eigenfrequencies computed (as in RefP^J) for 
our experimental aspect ratio L/H ~ 1.22. This shows 
that the translation of the grid, because of its specific 
drag profile, selects only a small set of specific inertial 
modes among the dense spectrum of modes evidenced 
in RefP^ In the next subsection, the spatial structures 
of the observed modes will be confirmed to match the 
theoretical structures predicted for the [n, to, s] modes 
of Table [I] Those selected modes are among the lowest 
order modes (n — 1 to 3), probably because their sim- 
pler structure better matches the spatial features of the 
ensemble-averaged flow at early times, i.e. before it is af- 
fected by rotation, and also because higher order modes 
decay faster by viscous damping. 

In order to check whether the inertial modes are 
present only in the ensemble- averaged flow, we also in- 
troduce the temporal Fourier transform of the total, 
u. a (x,z), and of the turbulent velocity field, xi! a (x, z), 
similarly to Eq. Q, from which we can define the corre- 
sponding energy spectra: 



Eturb{<y) 



1 

2^ 
1 

2V 



(8) 
(9) 



As for the kinetic energies, the energy spectra are addi- 
tive, E tot (a) = E. mean (a) + E turb (a). For a given total 
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FIG. 4. (Color online) Temporal energy spectrum of the total 
(dashed), mean (continuous) and turbulent (dashed-dotted) 
component of the flow as a function of <j/2f2 computed from 
40 decay realizations performed at Q = 0.84 rad s , with (a) 
linear and (b) logarithmic y-axis. Inertial modes can develop 
for pulsations a < 20.. The modes corresponding to the peak 
frequencies are given in Table [I] 



mode a/2Q 



[n, m, s] num. exp. 
[1,4,+] 0.2992 0.29 
[1,1,+] 0.4890 0.49 
[1,1,-] 0.6328 0.64 
[2,1,+] 0.7429 0.74 
[3,1,+] 0.8557 0.86 



TABLE I. Numerical values of normalized frequencies a/2Q 
for different modes of order [n, m, s] (where n is the vertical 
wavenumber, m characterizes the horizontal structure, and s 
is the symmetry of the mode) compared to the experimental 
peaks in Fig. [4] The uncertainty of the experimental values 
is ±0.01. The numerical values are computed for an aspect 
ratio identical the experimental one L/H = 1.22. 



spectrum E tot (a), the two spectra E mean and E tur b allow 
us to distinguish between two kinds of inertial modes: (i) 
reproducible (i.e. phase coherent) modes, characterized 
by peaks in E mean but not in E tur b\ (ii) non-reproducible 
(i.e. with random phase) modes, characterized by peaks 
in E tU rb but not in E mean . Interestingly, Fig. [4] shows 



the complete absence of peaks in the turbulent spectrum 
Eturb(&)i indicating that all the inertial modes present in 
our system are reproducible, and therefore not turbulent. 
This confirms that inertial modes in this system are more 
likely excited by reproducible flow features induced by 
the grid translation. This also suggests that the initial 
turbulence right after the grid translation is effectively 
not affected by the rotation, as could be expected from 
its relatively large grid Rossby number. Therefore, this 
initial turbulence is expected to be similar to classical 
grid turbulence in non-rotating systems. 

Finally, it is important to note that more than 97% of 
the energy of the mean flow lies in the range of pulsations 
[0, 20], confirming that almost all its energy is stored in 
inertial modes. This result suggests that the ensemble 
averaged flow is correctly converged, so that the spec- 
trum of the mean flow is almost not contaminated by the 
spectrum of residual turbulent fluctuations which spans 
over larger pulsations. On the contrary, we see that there 
is a significant amount of energy for a > 20 in the turbu- 
lent spectrum, which actually corresponds to the rapid 
small scales of usual 3D turbulence which are not directly 
affected by the rotation. 

C. Spatial structure of the inertial modes 

By performing a band-pass filtering of the mean ve- 
locity field at the frequencies of the peaks identified in 
Emean it is possible to extract the spatial structure 
of each inertial mode. However, this band-pass filter- 
ing must take into account the decay of the amplitude 
(inms(t) of each mode, which is actually reflected in the 
width of the peaks. Inserting Eq. ^ into Eq. (JoJ) yields 

U CT (x)= r"«™ s Wv„™ s We ,(u """- ff)t df. 

n,m,S Twin 

( 10 ) 

If the selected frequency a coincides with one of the peak 
frequency io nms identified in the spectrum, the extracted 
field is then equal to the spatial structure of the corre- 
sponding mode [n,m,s] weighted by its average ampli- 
tude, 

XJ Unms (x,z) = v nm8 (x,z) / a nms (t)dt. (11) 

JT min 

The resulting spatial structure U Wnms (x, z) is a com- 
plex vector field, whose real and imaginary parts encode 
the geometric features (amplitude, ellipticity, orientation 
and phase) of the elliptic orbit described by the velocity 
vector at each location during one period. Those or- 
bits can be simply reconstructed by plotting the trace 
of the real oscillating vector field, 3?(U Wnms e lLJ " mst ), for 
t £ [0, 2iv/u) nms ]. The resulting plots for the 5 dominat- 
ing inertial modes listed in Table [I] are shown in Fig. [5] 
The velocity field, taken for an arbitrary phase of the 
oscillation, is also illustrated by the vector arrows. The 
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color of the ellipses corresponds to the local ellipticity 
— 1 < e < 1, where positive ellipticity corresponds to an 
anticlockwise rotation. 

The geometrical features of these fields provide a good 
qualitative support for the identification of the frequen- 
cies given in Table [I] In particular, the modes of larger 
energy are of vertical wavenumber n — 1 [Figs. [5ja-c)], 
i.e. they show one cell in the vertical direction. Higher 
order vertical modes, n — 2 and n — 3, are also found, 
with much weaker energy. The irregular shape of the 
last mode, n = 3, is probably due to its very low en- 
ergy level, and to a possible mixing with other modes of 
similar frequency. Except for the antisymmetric mode 
[1, 1, -] at a/2tt = 0.64 [Fig.^c)], all the excited modes 
are symmetric (s = +) with respect to the rotation axis. 
This weak antisymmetric mode is probably excited from 
a spontaneous symmetry breaking of the flow, or from a 
slight residual left-right asymmetry of the grid. 

The horizontal structure of the modes cannot be in- 
ferred from these fields in the vertical plane. However, 
the visualizations provided in Ref™ indicate that, in 
spite of the square shape of the container, the modes 
are approximately axisymmetric, at least at a moderate 
distance from the rotation axis. 

The selection of the excited inertial modes among all 
the possible resonant ones may originate from a residual 
inhomogeneity of the jet velocities through each mesh of 
the grid. Such inhomogeneity is not necessarily associ- 
ated to a defect in the grid, but could arise from the influ- 
ence of the side walls of the tank or from a spontaneous 
symmetry breaking of the flow. In particular, the spatial 
structure of the two dominant modes (at a/2£l = 0.29 
and 0.49) suggests that larger (or smaller) velocities are 
produced near the center of the grid, initiating a marked 
vertical columnar oscillation. 

In order to check the possible inhomogeneity of the 
flow produced by the grid, we have measured the ve- 
locity profile just behind the grid during its translation 
(Fig. [6^,, and the continuous curve in Fig. [6b). Slightly 
weaker jets are indeed encountered near the borders of 
the grid, probably originating from the friction with the 
sidewalls. The resulting rounded average profile could be 
responsible for the excitation of the dominant symmetric 
modes. 



IV. MODIFIED CONFIGURATION WITH 
INNER WALLS 

A. Comparison with a previous configuration 



FIG. 5. (Color) Spatial structure of the 5 dominating inertial 
modes listed in Table [I] extracted by band-pass filtering of 
the ensemble-averaged fields. The ellipses show the velocity 
orbit, and the arrows illustrate the velocity field at a given 
arbitrary phase of the oscillation. The color of the ellipses 
traces the ellipticity — 1 < e < 1 (see the text for details), 
(a) [1,4,+]; (b) [1,1,+]; (c) [1,1,-]; (d) [2,1,+]; (e) [3,1,+]. 
Resolution of the fields has been reduced by a factor 6 for a 
better visibility. 



Although the excitation mechanism of inertial modes 
in a rotating grid turbulence experiment is interesting 
in itself, generating such turbulence without those re- 
producible modes in a confined geometry is also highly 
desirable. This is motivated by the fact that most the- 
oretical and numerical works on this subject have been 
performed in purely homogeneous rotating turbulence, 
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FIG. 6. (Color) Vertical velocity field in the wake of the grid, 
(a) Simple grid configuration, (b) Modified configuration with 
the inner walls attached to the grid, represented by the 6 ver- 
tical thick lines. The grid is towed from the bottom, and is 
at the height z = 290 mm in these snapshots, (c) Horizontal 
profile of the vertical velocity, at a distance 8z — 5 cm below 
the grid (indicated by the horizontal lines in a and b). Con- 
tinuous line: simple grid configuration. Dashed line: modified 
configuration with inner walls. The vertical black lines show 
the locations of the vertical inner walls. 



without inertial modes (see, e.g., Sagaut and CamborPSI 
and references therein). In this section, we character- 
ize the flow generated by a modified grid configuration, 
which turns out to be almost free of inertial modes. 

In a similar rotating turbulence experiment, in order to 
avoid a large scale recirculation flow triggered by the grid 
translation, Staplehurst et affl have designed an original 
setup in which a parallelepipedic inner tank, consisting 
in 4 vertical sidewalls without top and bottom walls, was 
attached on the upper side of the grid. This inner tank, 
of 35 cm side, was translated with the grid from the top 
to the bottom, inside the primary 45 cm wide square con- 
tainer filled with water. In this situation, the inner side- 
walls are therefore located downstream the grid, and once 
the grid is lowered, the effective flow volume is reduced 
to the size of the inner tank. With this additional in- 
ner tank, these authors actually succeeded to reduce the 
recirculation flow by a factor 5. However, only the non- 
rotating case was described, and the possible presence of 
inertial modes with rotation has not been investigated. 




FIG. 7. (Color online) Schematic view of the modified grid 
configuration (the outer water tank is not shown; see Fig. [TJ . 
Three inner tanks are mounted on the grid, and turbulence 
is generated by raising the set grid+tanks. Each inner tank 
consists in 4 vertical sidewalls, without top and bottom walls. 



The purpose of the inner tank in the previous down- 
stream configuration was to block the transverse motions 
in the wake of the grid. Here, in order to force equally 
distributed jets through the grid, while keeping the full 
volume for the working fluid, we propose to modify the 
solution of Staplehurst et aZ.pby adding inner walls lo- 
cated upstream the grid. In our configuration, the inner 
walls are still attached on the upper side of the grid, but 
the grid here is lifted from the bottom to the top. As 
a consequence, the working fluid volume in the wake of 
the grid is the same as in the simple grid configuration, 
so that direct comparisons of the flow structure can be 
performed between the two configurations. Several con- 
figurations have been tested, with 1, 2 and 3 inner tanks. 
Best results have been obtained with 3 inner tanks, and 
only this configuration is described in the following. 



B. Modified experimental setup 

The new setup is shown in Fig. [7j It consists in three 
parallelepipedic PVC tanks, of respective height 50, 40 
and 40 cm and width 40, 24 and 8 cm. Each tank con- 
sists in 4 vertical sidewalls, without top and bottom walls. 
The 3 tanks are fit together co-axially, so that two neigh- 
bour sidewalls are separated by two grid meshes, and 
attached on the upper side of the grid. All other param- 
eters, in particular the grid velocity and Rossby number, 
are the same as for the simple grid configuration. Since 
our grid is translated from the bottom to the top, here 
the inner tanks are upstream the grid, so that the decay- 
ing flow is not confined inside the smallest tank. 

The presence of the inner tanks implies that no hori- 
zontal cover could be placed to define the upper bound- 
ary of the flow, contrarily to the "simple" grid configu- 
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FIG. 8. (Color online) Total (dashed), mean (continuous) 
and turbulent (dashed-dotted) kinetic energies as a function 
of reduced time Qt/2-K from 40 decay realizations performed 
at Q. = 0.84 rad s _1 with the modified configuration with 
inner tanks. Inset: ratio of turbulent to total kinetic energy 
as a function of reduced time Qt/2ir with (points) and without 
(triangles) inner tanks. 



ration. However, we can consider here that an effective 
upper boundary is approximately defined by the grid it- 
self, locked at a height of 46 cm after its translation. In- 
deed, we have checked that the "simple" grid configura- 
tion without upper rigid boundary also produces inertial 
modes, which are as intense as with the rigid boundary. 
Moreover, here, the inner walls efficiently block the flow 
in the region between the grid and the free surface. 

The velocity field and profiles in this modified configu- 
ration are shown in Figs. [6](b) and (c), compared to that 
obtained with the simple grid configuration [Figs. [6^a) 
and (c)]. The main difference is that, now, a stronger 
upper flow is present in the wakes of those grid bars sup- 
porting the sidewalls, whereas the upper flow in the wake 
of the other grid bars is significantly reduced, and actu- 
ally even difficult to detect with our PIV resolution. As a 
consequence, two neighbor jets originating from the same 
"passage" between two sidewalls tend to gradually merge. 
At this point, comparing the velocity profiles of the two 
configurations, with and without the inner tanks, it is 
not clear whether changes in the energy budget between 
turbulence and mean flow may arise from the observed 
differences. 



C. Turbulence generated with the grid + inner 
tanks configuration 

The decay of the kinetic energies of the total, mean 
and turbulent components of the flow as a function of 
time, from ensemble averages over 40 decay realizations 
with the modified configuration, is shown in Fig. [8] In 
comparison with the simple grid configuration (Fig. [3]), 
the oscillations due to the inertial modes are strongly re- 
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FIG. 9. (Color online) Temporal energy spectrum of the total 
(dashed), mean (continuous) and turbulent (dashed-dotted) 
component of the flow as a function of a/20, computed from 40 
decay realizations performed at £1 = 0.84 rad s with inner 
tanks, with (a) linear and (b) logarithmic j/-axis. The dotted 
line reproduces spectrum of the mean flow for comparison 
(Fig. HJ, obtained with the simple grid configuration. 



duced, both for the mean and the total kinetic energies. 
We also see that, contrary to the simple grid configu- 
ration, the turbulent kinetic energy is now significantly 
larger than the one of the mean flow. This can be bet- 
ter seen in the inset of Fig. |8j showing that turbulence 
contains, after a transient of about one tank rotation, ap- 
proximately 85 ± 5% of the total kinetic energy, a value 
much larger than the 50% obtained with the simple grid 
configuration. 

In Fig. [9j we superimpose the temporal energy spec- 
trum of the mean flows obtained with (continuous line) 
and without (dotted line) the inner tanks. We clearly see 
that, with the inner tanks, the two dominant peaks (at 
a/2fl — 0.29 and 0.49) have been reduced by more than a 
factor of 10. The reduction of the other peaks, in partic- 
ular the n — 2 mode at <j/2VL ~ 0.74, is less pronounced. 
The decrease of the two dominant inertial modes essen- 
tially explains the significant drop of the kinetic energy 
of the mean flow with respect to turbulence shown in 
Fig. [8] We can also note slight frequency shifts between 
the two configurations, probably originating from the not 
perfectly identical boundary conditions at the top of the 
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container, and hence from slightly different effective as- 
pect ratios. 

Interestingly, incrtial modes are not only absent from 
the mean flow spectrum, E mean {cr) but also from the 
turbulent spectrum E tU rbis T ). This indicates that their 
disappearance in the modified grid configuration cannot 
be attributed to a change in the nature of the excited 
inertial modes, from reproducible to non reproducible, 
but to a true disappearance. 

The differences in the energy decays and spectra show 
that the addition of the inner tanks to the grid leads to 
a much more efficient transfer of the initial mean flow 
energy into turbulence, bypassing the generation of iner- 
tial modes. However, the origin of the changes between 
the two configurations is not clear. In particular, it is 
not evident if, and how, the differences in the jet profiles 
in Fig. [6^c) may explain the preferential energy transfer 
from the mean flow to the turbulence in the modified con- 
figuration. These results illustrate the sensitivity of the 
inertial modes production to slight changes in the grid ge- 
ometry, at least for the grid Rossby number considered 
here, Ro g = 10.4. 

V. DISCUSSION AND CONCLUSION 

In this article, we characterize in details the flow gener- 
ated when a simple grid is rapidly translated in a rotating 
and confined volume of fluid. We show that, in addition 
to a non-reproducible turbulent flow, the grid translation 
initiates a reproducible mean flow composed of resonant 
inertial modes, which contains a significant amount of the 
total kinetic energy. This observation agrees with the re- 
cent experiments of Bewley et alW in a similar geometry. 
In addition to quantitative comparisons of the mode fre- 
quencies with the numerical predictions of Maas^l, we 
also provide here good qualitative comparisons for the 
spatial structure of these modes. 

These results suggest that the turbulent (i.e. the non 
reproducible) component of the flow generated in this 
configuration cannot be considered as freely decaying. In- 
deed, energy transfers between the inertial modes and the 
turbulence may exist, so that the energy initially stored 
in the inertial modes may continuously feed the turbu- 
lence. However, measuring this transfer would require 
to compute the turbulent Reynolds stress tensor and its 
coupling with the ensemble-averaged flow. This cannot 
be performed from the PIV data presented here, which 
are limited to two velocity components in a single plane. 

The ability of a translated grid to generate or not re- 
producible inertial waves or modes depends on the grid 
Rossby number, Ro g = V g /2flM, and also on the geo- 
metrical details of the grid configuration, in particular 



its solidity (ratio of solid to total area). Indeed, for a 
given grid geometry, a low Rossby number seems more 
likely to excite inertial modes, because the reproducible 
wake pattern of the grid, or the induced reproducible re- 
circulation flow, may directly force inertial modes with a 
phase coherence. On the other hand, we could speculate 
that, for much larger Rossby numbers (and assuming that 
the Reynolds number is very large too), the reproducible 
wake pattern has more time to loose its coherence before 
being affected by rotation, so its energy may be more ef- 
ficiently transferred to turbulence than reproducible in- 
ertial modes. Even if inertial waves are generated from 
these decaying incoherent motions, their random phase 
should cancel them out in the ensemble average, so iner- 
tial modes should not be found in the mean flow. 

In the experiments of Bewley et aL^the moderate grid 
Rossby number Ro g of 5.5 is found to produce strong in- 
ertial oscillations. In the experiments of Staplehurst et 
aZ.p Ro g lies in the range 1 — 3.2, so that significant 
inertial modes are also expected, although they are not 
described by these authors. Finally, in the present exper- 
iment, for a slightly larger grid Rossby number of 10.4, 
the production of reproducible inertial modes is found 
to be very sensitive to slight details in the geometry of 
the grid. Indeed, by adding a set of inner tanks to the 
grid, we show that it is possible to significantly reduce 
the production of incrtial modes, although the mecha- 
nism responsible for this preferred energy transfer to- 
wards turbulence instead of modes is not elucidated in 
the present work. This sensitivity of the flow evidenced 
for this particular Ro g is consistent with the results of 
Morize and MoisyP^ covering a range of Ro g between 2 
and 65, in which a smooth transition around Ro g ~ 10 
is found in the decay law of the turbulent energy, sug- 
gesting a change in the energy transfer between inertial 
modes and turbulence. 

From those results, one may conclude that, providing 
the grid Rossby number and the geometrical features of 
the grid are carefully selected, achieving a nearly "pure" 
rotating turbulence, free of reproducible inertial modes, 
is actually possible in a confined geometry. This system is 
therefore suitable to explore experimentally the influence 
of the rotation on freely decaying turbulence. 
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